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Abstract 

It is shown that the coupling of heat with acceleration first proposed by Eckart would have 
an overwhelming effect in the growth of density mass fluctuations, even in non-relativistic 
fluids in the presence of a gravitational field. Gravitational effects would be negligible if the 
heat-acceleration relation is assumed to be valid for the hydrodynamic equations. A direct 
implication of this result is that recent alternative first order in the gradients theories must 
be taken into account while describing a special relativistic fluid. 
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I. INTRODUCTION 



Astronomical observations strongly suggest that structures are continuously being 
formed in our Universe. According to our current knowledge, gravity is the leading 
interaction responsible for structure formation. The critical parameters in these 
events are the particle density and certain lengths that have been known since the 
pioneering works of J. Jeans |l| U Also, the characteristic formation time that can 
be inferred from these theories lead to very large values compared with those that 
appear in microscopic processes such as particle collisions. 

In the original works regarding structure formation, the growth of linear density 
perturbations has been analytically modeled in terms of the basic conservation laws 
of particle number and momentum. In this approach, density number perturbations 
satisfy a wave equation whose possible solutions provide basic values for the critical 
parameters that allow density growth, or simply acoustic wave propagation |2|. On 
the other hand, dissipative processes are always present in real phenomena. The effect 
of viscosity and heat flux on the determination of the criticalparameters present in 
gravitational collapse have been addressed by several authors p| J4| most of those 
works use a classical (non relativistic) framework in order to include dissipation. 

In the present work we show that the inclusion of the heat-acceleration coupling 
in the framework of relativistic hydrodynamics would yield a pathological behav- 
ior regarding the observed characteristic times in which gravitational collapse takes 
place. Indeed, the generic instability first found by Hiscock and Linblom |6|] and later 
thoroughly examined in a physical context [7| would have a dominating (unobserved) 
effect over gravity. The paper is divided as follows, section two includes the basic 
equations in which Eckart's framework is developed, including weak gravitational ef- 
fects; in section three the evolution of the first order fluctuations that arises from the 
linearized equations is analyzed. The corresponding dispersion relation is examined 
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in the presence of gravity in order to establish the dominant effects. Finally, section 
four is devoted to a discussion of the results here obtained and its corresponding 
implications. 



II. RELATIVISTIC TRANSPORT EQUATIONS IN ECKART'S FRAME 

The starting point in Eckart's formalism for a simple fluid is the particle conser- 
vation equation 

% = 0. (1) 

In Eq.([IT), for a number density n and hydrodynamic velocity w M , = nu^ denotes 
the particle flux. In terms of total derivatives, this equation reads 

h + n0 = (2) 

where 9 = uf" . The energy- momentum tensor balance has the form 

= (3) 

where 

TIF I I 

T* = —u»u v + ph$ + + -q»u v + -u»q v (4) 

c c c 

Here, e is the internal energy per particle, p is the local pressure, and h% = 5% + 

is the spatial projector with signature (+ + +— ). The Navier tensor IT^ and the heat 

flux g M satisfy the orthogonality relations: 

u^W v = u v W v = 0, ?X = 9% = ( 5 ) 

Now, the projection u u T£^ = leads to the internal energy balance equation, 

ne + p9 + + + \u v q v = (6) 
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and, together with the expression T^.^ = 0, to the equation of motion: 

fne p \ . fne p \ , , . _.. 

+ ^ (q^u u + q^u^ + 9q u + u"^) = (7) 

Eq. ([6]) can be rewritten in terms of n and T using the functional hypothesis e = 
(In) T ™ + \w) T l ea dirig to the expression 

nC n f + p6 + u^W u + g£ + -^g„ = (8) 

The entropy balance is also established by means of the functional hypothesis s = 
(in) ^ + (If) l ea ding to the expression for the entropy production a [8]: 

This last equation has motivated the introduction of the heat-acceleration coupling 

= -Kh; (t* + (10) 

where K is the thermal conductivity coefficient. In the next section we shall analyze 
critically the implications of Eq. ffTOl) using the linearized version of the set of 
equations ( 121171) and Eq. ([8]) in the context of special relativity, and introducing a 
gravitational field. 



III. LINEARIZED EQUATIONS: EXPONENTIALLY GROWING MODES 

In order to examine the behavior of small fluctuations 8X around the equilibrium 
values of a thermodynamical variable X, we shall write X = X + 5X, where the 
subscript o denotes the average value. Also, in the right hand side of Eq. ([7J we will 
explicitly write a gravitational force, taking into account a Newtonian approximation 
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for the covariant derivatives included in this expression. Neglecting second order 
terms in the fluctuations, the resulting equations in the comoving frame u = read: 

^(5n)+n 59 = (11) 

1 c) 

-a (n e + Po ) —(86) + kT a V 2 {5n) + n W 2 (5T) - AV 2 {59) 

Cr Ot 

-4v4(*o - ^^.m = i + po) v 2 (m (12) 



nCJT + g?) ^ - k (<5T«* + ^<5u fc ) = 0. (13) 



Use has been made of the fact that p = p(n, T). In Eq. (112"]) we have introduced 
the gravitational potential ip, which satisfies in the weak field limit the Poisson equa- 
tion V 2 5(p = —A7rGm5n. Also, the linearized expression of Eq. (J2J) has been written 
in terms of the longitudinal mode 69. 

In the Fourier-Laplace space, the system of Eqs. (I11H13P reads: 

s5n + n o 30 = 5h(q, 0) (14) 

1 ~ ~ kT ~ 

- q 2 kT 3h + ATcGm 2 n 3h + — (n e + p a ) s39 + Aq 2 39 ^s 2 39 

& c 4 

-q 2 n k3T + ^-s3T = ~^s39(q, 0) + \ (n e o + Po ) 59(q, 0) - ^ST(q, 0) (15) 

c c c c 

~ kT ~ ~ ~ ~ D uT 

kn T S9 £s50 + c n n sST + nq 2 8T = ST(q, 0) + -^^59(q, 0) (16) 

c c 



In Eqs. (I14H IB"]) , for a given thermodynamical fluctuation SX(f, t), the correspond- 
ing Fourier-Laplace transform is denoted as SX = SX(q, s). 
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The dispersion relation associated to the system (I14H16P governs the dynamic 
behavior of the thermodynamic fluctuations, and read: 



s n 

-kT Q q 2 + 4nGm 2 n p Q s + Aq 2 + ^s 2 -q 2 n Q k + ^-s 

kn T - +c n n s + nq 2 

The resulting dispersion relation has the form: 

a s 4 + a t s 3 + a 2 s 2 + a 3 s + a 4 = 



where 



«o 



PoC 4 



;i7) 



;is) 



(19) 



2 2kn q 2 T n 2 _ 
a 2 = -4c„n g h <? «p 



(20) 
(21) 



a 3 = -47iGM 2 c n nl + c n kn 2 T q 2 + k 2 n 2 T q 2 + Anq 



2rp „2 , /„2„2t-t „2 



(22) 



a 4 = — 47rGM n G g /c + kn T q k 



(23) 



The gravitational effect only appears in Eqs. (22-23). In this case, the denomina- 
tor c 4 in Eq. (19) generates a very large root of the full dispersion relation that can 
be accurately approximated by solving the simple equation: 

a Q s A + ais 3 = (24) 

The non-trivial root of Eq. 024p is the one obtained in reference [6| . It is relevant to 
notice that, in the absence of dissipation k — 0, A=0, the dispersion relation reduces 
to a simple cubic equation that leads to the well-known Jeans wave number. 
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The conclusion of this analysis is that the heat acceleration coupling would lead to 
extreme short time scales (of structure formation even in the presence of gravitation 



10 34 seconds for water at 300 K and one atmosphere of pressure j(3|). This is a 
serious drawback for Eckart's coupling Eq. ( TTO]) . 



IV. FINAL REMARKS 



This work supports the points raised in earlier work 12|] where it was shown that 
the use of the heat-acceleration coupling, introduced since 1940, yields unobserved 
predictions regarding generic instabilities and light scattering spectra. The results 
here obtained show that this kind of coupling would also generate a pathological 
behavior regarding structure formation. Indeed, the introduction of Eq. (fTUj) in the 
balance equations leads to negligible effect of gravity. This is of course in contra- 
diction with current knowledge. Moreover, the introduction of the so-called second 
order theories may modify significantly the results here. In most version of those 
theories the corresponding constitutive equation reads, neglecting viscous effects: 

^ + T ? = - Kh ; (V + (25) 

where r is relaxation time much larger the characteristic time of 10~ 34 seconds. In 

this equation the heat acceleration coupling is still retained. On the other hand, a 
possible first order alternative constitutive equation has been successfully introduced 
in Ref. [?']. The study of these alternatives regarding the gravitational instability in 
the presence of these constitutive equations will be addressed in the near future. 
The authors acknowledge D. Brun-Battistini for her valuable comments for this 
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